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Abstract 



We introduce a new graph polynomial that encodes interesting properties of graphs, for 
example, the number of matchings and the number of perfect matchings. Most importantly, 
for bipartite graphs the polynomial encodes the number of independent sets (#BIS). 

We analyze the complexity of exact evaluation of the polynomial at rational points and 
show that for most points exact evaluation is #P-hard (assuming the generalized Riemann 
hypothesis) and for the rest of the points exact evaluation is trivial. 

We conjecture that a natural Markov chain can be used to approximately evaluate the 
polynomial for a range of parameters. The conjecture, if true, would imply an approximate 
counting algorithm for #BIS, a problem shown, by [TO], to be complete (with respect to, so 
called, AP-reductions) for a rich logically defined sub-class of #P. We give a mild support for 
our conjecture by proving that the Markov chain is rapidly mixing on trees. As a by-product we 
show that the "single bond flip" Markov chain for the random cluster model is rapidly mixing 
on constant tree-width graphs. 

1 Introduction 

Graph polynomials are a well-developed area useful for analyzing properties of graphs (see, e.g., 
the following survey papers [131 114| ). Arguably the most intriguing graph polynomial is the Tutte 
polynomial 38 , 39 . The partition function of the random cluster model from statistical mechanics 
provides a particularly simple definition: for a graph G — (V, E) let 



where k(S) is the number of connected components of the graph (V,S). (It is well-known that 
the Tutte polynomial is obtained from Z by a simple transformation, see, e.g., equation (I12[) in 
Section [3]) The Tutte polynomial includes many graph polynomials as special cases, e.g., the 
chromatic polynomial, the flow polynomial, and the Potts model (see, e.g., [H]). 
Now we define our graph polynomial. 

Definition 1. The R2- polynomial of a graph G — (V,E) is 
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where rk^S) is the rank of the adjacency matrix of (V, S) over ¥ 2 (the field with 2 elements). 

Now we look at how R 2 (G; q, //) encodes some properties of graphs. For q — pr 1 ! 2 equation ^ 
becomes 

P(G;ri:=R 2 (G;n- 1/2 ,ri= £ ^ k ^. 

SCE(G) 

We claim that P(G; 0) is the number of matchings in G. To see this, note that rk 2 (S) < 2\S\ (since 
adding an edge to S changes two entries in the adjacency matrix and hence can change rank by at 
most two), and ik 2 (S) < 2|5| if S is not a matching (since rank of the adjacency matrix of a star 
is 2 < 2 1 5 1, and adding further edges preserves the strict inequality). 
Let 

P(G;t,/ J L):=t^R 2 (G;l/t,fj,) and P 2 (G; /1) := ^ |y|/2 P(G; 0, fj,). 

Then P 2 (G;0) is the number of perfect matchings of G. To see this, note that that only subsets 
with full rank adjacency matrix contribute to P(G;0,/i), and then only the minimal cardinality 
subsets with full rank adjacency matrix contribute to P 2 (G; 0) (these subsets are exactly the perfect 
matchings) . 

From now on we focus solely on bipartite graphs. For a bipartite graph G = (U U W, E) we let 

R' 2 (G;X,fi)= ]T A rka < s V |S| , (3) 

SC_B 

where rk2(£) is the rank of the bipartite adjacency matrix of (U U W, S). Note that 

i? 2 (G;A, M ) =R' 2 (G;\ 2 ,v), (4) 

since the adjacency matrix contains "two copies" of the bipartite adjacency matrix (one of them 
transposed). (The reason for definition ([3]) is that we prefer to operate with bipartite adjacency 
matrix for bipartite graphs.) 

In Section [5] we prove that R' 2 counts the number of independent sets in bipartite graphs. 

Theorem 1. Let G = (U U W,E) be a bipartite graph. The number of independent sets of G is 
given by 

2 \u\+\w\-\e\ r ^ g . 

We show in Section[3]that exact evaluation of the polynomial R' 2 (G; A, /i) is #P-hard at a variety 
of rational points (A,/z) assuming the validity of the generalized Riemann hypothesis (GRH). 

Theorem 2. Exact evaluation of R' 2 at rational point (A, (i) is 

• polynomial-time computable when A £ {0, 1} or fi ~ or (A,/i) = (1/2, — 1); 

• #P-hard when A ^ {0, 1, 1/2} and /j ^ 0, assuming GRH; 

• #P-hard when A = 1/2 and /x ^ {0, —1}. 

Remark 1. For the non-bipartite case we have the following classification. Exact evaluation of R 2 
at rational point (X, n) is polynomial-time computable when \i — or A G { — 1,0,1}; the A = — 1 
case follows from the fact that a skew- symmetric matrix with zero diagonal has even rank over 
any field (the zero diagonal condition is redundant for fields of characteristic ^2). For any other 
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rational X and /i we get #P-hardness of evaluating the R2 polynomial from Theorem [D and ([¥]) 
(again assuming GRH). (Note that (A,/i) i-> (A 2 ,/i) never maps to the easy case (1/2,-1), since 
X is rational. It would be nice to have hardness classification of evaluating R2 and R' 2 for, say, 
algebraic A and jj,.) 

Because of the hardness of exact evaluation of R' 2 , we turn to approximate evaluation of 
R' 2 (G;X,»). 

Remark 2. There is a fully polynomial randomized approximation scheme (FPRAS) for R' 2 (G; A,/i) 
when A = 1/2 and —1 < fj, < 0. This follows from Theorem\5\ (stated in Section^ and the fact 
that there is an FPRAS for #PJ5IS(rj) (defined in Section^ when < 77 < 1 (see [19], and note 
that #PBIS(?7) corresponds to f3 = 1, 7 = (1 + rj)/(l — rj) and in their parametrization) . 

We now define the sampling problem associated with R' 2 . 

Rank Weighed Subgraphs with \,fi > 0, (RWS(A,/^)) 
Instance: a bipartite graph G = (U U W, E), 
Output: S C E with probability of S cx A rk2 ( s V' 5 '- 

The "single bond flip" chain is a natural approach to sampling from RWS(A,/^). 

Definition 2. Single bond flip chain is defined as follows: pick an edge e £ E at random and let 
S = X t ® {e}. Set X t+ i = S with probability (1/2) min{l, \^2(s)-rk 2 (x t ) ^\s\-\x t \} md Xf+i _ x± 
with the remaining probability. 

In each step of the single bond flip chain, we have to compute the rank of a matrix over F2 
(corresponding to S) which differs from the current matrix (corresponding to X t ) in a single entry. 
One can use dynamic matrix rank problem algorithms to perform this computation in 0(n 1575 ) 
arithmetic operations per step [16] . 

Instead of flipping one edge in a step, we can have another Markov chain which flips a random 
subset of edges adjacent to a single vertex. It seems likely that the new chain can generate good 
random samples faster than the single bond flip chain — a step of the new chain can be performed in 
0{n 2 ) arithmetic operations (using "rank one update" for the dynamic matrix rank problem [16]). 

We optimistically conjecture that the single bond flip chain (and hence the other chain) mixes 
rapidly. 

Conjecture 1. For any fixed X, fi > the single bond flip chain mixes in polynomial time. 

Note that Goldberg and Jerrum [T5] conjecture the opposite, in particular for A = 1/2 and 
/i = 1. We make the conjecture based on similarity of R' 2 polynomial to the random cluster model 
and the fact that slow mixing of the single bond flip chain was not established in random cluster 
model (slow mixing is usually easier to establish than rapid mixing). 

Conjecture Q] could be false for small A and large pL if the problem of maximizing | S'l — c • rk 2 (5') 
(c is a constant and c > 1) over S C E is hard (the corresponding maximization problem is easy 
for the random cluster model). 

We can provide a mild support for Conjecture [T] In Section [S] we prove that for fixed A, [i > 
the single bond flip chain mixes, in time polynomial in the number of vertices, for trees. 

Theorem 3. For every fixed A, fi > 0, the mixing time r(e) of the single bond flip chain for a tree 
on n vertices is 

t{e) = O (n 3+ l lo ^ A l(|logA| + |log/i| +log(l/e))) . 
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As a by-product of our techniques, we show that single bond flip Markov chain for the random 
cluster model is rapidly mixing if q, fi > and G has constant tree-width (the condition q, fi > is 
equivalent to x,y > 1 for the Tutte polynomial T(G; x,y)). 

2 Independent sets in bipartite graphs 

The problem of counting independent sets (#IS) in a graph is of interest in both computer sci- 
ence and statistical physics (independent sets are a special case of the so-called hard-core model, 
see, e.g., [S]). Exact computation of #IS is #P-complctc even for 3-regular planar bipartite 
graphs 40, 45 . Fully polynomial randomized approximation scheme (FPRAS) is known for graphs 
with maximum degree A < 5, [28J [T2l [43] . Unless RP=NP, an FPRAS does not exist for graphs 
with A > 25, [5] . 

Now we focus on the problem of counting independent sets in bipartite graphs (#BIS). While 
for exact counting the complexity of #BIS and #IS is the same, the situation looks very different 
for approximate counting, for example, no inapproximability result is known for #BIS. Dyer et 
al. [10] show that #BIS is complete w.r.t. approximation-preserving reductions (AP-reductions) in 
a sub-class of #P. Many problems were shown to be equivalent (w.r.t. AP-reductions) to #BIS, 
for example, #Downsets, #1p1nSat [ID] , computing the partition function of a ferromagnetic 
Ising model with local fields [17] , and counting the number of satisfying assignments of a class of 
Boolean CSP instances [IT]. A pertinent negative result for #BIS is that Glauber dynamics (or 
more generally, any chain whose states are independent sets and that flips at most 0.35n vertices 
in one step) cannot be used to efficiently sample random independent sets in a random 6-regular 
bipartite graphs onn + n vertices [§] . 

The rest of this section is devoted to proving Theorem [T] It will be convenient to work with 
matrices instead of graphs. For two zero-one matrices A, B we say B < A if B corresponds to a 
subgraph of A, formally 

Definition 3. Let A, B be zero-one n\ x matrices. We say B < A if Aij = implies Bij = 0, 
for all i £ [ni] and j £ [712]- Let Ca be the set of zero- one n\ x n-i matrices B such that B < A. 

Let #i(A) denote the number of ones in A (that is, the number of edges in the corresponding 
graph). The RWS problem rephrased for matrices is: 

Rank Weighed Matrices with A,/i > (RWM(A,/x)) 
Instance: an n\ x 712 matrix A. 

Output: B (EC A with probability of B oc A rk2 ( B V #l(B) - 

The problem of sampling independent sets in bipartite graphs is: 

Bipartite Independent Sets (BIS) 

Instance: a bipartite graph G = (U UW,E). 
Output: a uniformly random independent set of G. 

Before we show a connection between BIS and RWM(l/2, 1) we remark that to sample bipartite 
independent sets it is enough to sample a subset of one side, say U, from the correct (marginal) 
distribution. We now describe this distribution in a setting which will be advantageous for the 
proof of Theorem [I] 

We will represent an independent set by a pair of (indicator) vectors u, v (where u € F^ 1 and 
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Definition 4. We say that two vectors a,/3 £ F£ share a one if there exists i £ [n] such that 
cm = ^ = 1. 

We will use the following simple fact. 

Observation 1. Let a, (3 £ ¥% ■ Let d be the number of ones in (3. 

• if a, (3 share a one then there are 2 d ~ 1 vectors f3' < f3 such that a T f3' = mod 2. 

• if a, (3 do not share a one then there are 2 d vectors f3' < (3 such that a T (3' = mod 2. 

Let u £ F'2 1 be a vector. We would like to count the number of v £ FJ? 2 such that u, v is an 
independent set. Note that u, v is an independent set iff Vj = for every j £ [112] such that u and 
j/'-th column of A share a one. Let k be the number of columns of A that do not share a one with u. 
Then we have 

u £ Fj 1 occurs in 2 k independent sets. (5) 

Thus to sample independent sets in a bipartite graph G with n\ x 712 bipartite adjacency matrix A 
it is enough to sample u £ V 1 ^ 1 with the probability of u proportional to 2 k , where k is the number 
of columns of A that do not share a one with u. We will call this distribution on u the marginal 
BIS distribution. 

The following lemma shows a tight connection between BIS and RWM(l/2, 1) — given a sample 
from one distribution it is trivial to obtain a sample from the other one. 

Lemma 1. Let G be a bipartite graph with bipartite adjacency matrix A. 

• Let u,v be a uniformly random independent set ofG. Let B be a uniformly random matrix from 
the following set {D £ Ca \ u T D = mod 2}. Then B is from the RWM(l/2, 1)- distribution. 

• Let B £ Ca be a random matrix from the RWM(l/2, 1) -distribution. Let u £ F^ 1 be a 
uniformly random vector from the left null space of B (that is, {(3 £ F^ 11 | (3 T B = mod 2}). 
Then u is from the marginal BIS distribution. 

Proof. Let Q be the set of u, B pairs such that u T B = mod 2 and B < A. Let i\) be the uniform 
distribution on Q. Note that ip marginalized over u yields the RWM(l/2, l)-distribution on B < A, 
here we are using the fact that a d-dimensional space (in this case the left null space of B) over F2 
has 2 d elements. Formally, 

1 2" 1 ~ 1 ' k2 (- B ) 2~ rk2 (- B ) 

P{B)= TR J ri2 M = ioi = b$g; iAi)' (6) 

u\u L B=0 mod 2 

Next we show that tp marginalized over B yields the marginal BIS distribution. We compute 
the number of B < A such that u T B = mod 2. Let us use the same k as in ([5]), that is, k is the 
number of columns of A that do not share a one with u. 

Note that the columns of B can be chosen independently and only if the column and u share 
a one is the number of choices (for that column) halved. Let #i(^4) be the number of ones in A. 
Thus 

there are 2* l(A ^ (n2 - k ) choices of B < A such that u T B = mod 2. (7) 
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Note that for fixed u the counts in ([5]) and (O differ by a factor of 2# 1 ^^ ™ 2 (which is independent 
of u). Thus ip marginalized over B yields the marginal BIS distribution on u. Formally 

2 #i(A)-(n 2 -fc) 2 k 

p(u) = — ioi — = #ms(G) • (8) 

Note that this proves both claims of the lemma since in both cases the it, B pair is from ip (by 
first sampling from a marginal and then sampling the remaining variable) and the conclusion in 
both claims is a statement about marginal (of the remaining variable). □ 

Theorem [TJ now follows from the proof of Lemma [TJ 

Proof of Theorem [7J Let Q be the set from the proof of Lemma [TJ From ([6]) we obtain 

\Q\=R' 2 (G; 1/2 A)2 n \ (9) 

From ([8]) we have that the number of independent sets of G is given by 

# BIS ( G ) = wSU- ( 10 ) 

Combining (j9j) and (|10|) we obtain the theorem. □ 

We do not know a good combinatorial interpretation for the mod-2 rank of B for general graphs. 
For forests (which are, of course, always bipartite) we have the following characterization. 

Lemma 2. Let G — (V,E) = (U U W, E) be a forest with bipartite adjacency matrix A. Then 
rk 2 (A) is the size of maximum matching in G. 

Proof. Let a £ V be a leaf of G and let e = {a, b} £ E be the edge adjacent to a. Note that b is 
matched in every maximum matching M (otherwise one could add e to M). Thus removing b and 
all adjacent edges decreases the size of maximum matching by 1. 

Now we argue that removing b (and all adjacent edges) also decreases rank (over F2) by 1. 
W.l.o.g. assume that b corresponds to the first row and a corresponds to the first column. Removing 
b (and all adjacent edges) corresponds to removing the first row of A. Note that this decreases rank 
by at most 1 and it does decrease it by 1, since the only non-zero entry in the first column is in the 
first row. □ 

3 Exact evaluation of R' 2 (proof of Theorem [2J) 

We will prove Theorem [5] in this section. Let G = (V, E) = (U U W, E) be a bipartite graph. First 
we deal with the cases where exact evaluation of R' 2 {G: A, fi) is easy. For cases A £ {0, 1} and fi = 
we have 

R! 2 (G;0,[ J ,) = R' 2 (G;\,0) = l and R' 2 {G- 1, /*) = (1 + 
For A = 1/2 and /j = -1 we will show 

R' 2 (G; 1/2,-1) =2^1-^1+*, (11) 

where t is the number of isolated vertices in G (see the remark after Theorem [5] in this section). 
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For the hardness results we give reductions from the problem of evaluating the Tutte polynomial 
(to establish the second part of Theorem [2]) and #BIS (to establish the third part of Theorem [2]). 
The Tutte polynomial of a graph G = (V, E) is a polynomial in two variables x, y defined by 

T(G;x,y) = J>- l^W-Ofo - i)I*I-WmW 

SCE 

where k(S) is the number of connected components of the graph (V,S). The Tutte polynomial is 
closely related to the random cluster model (see, e.g., jM]). Let Z(G;q,/j.) be defined as in (JTJ. 
One ticts 

T(G; x,y) = (x - l)-"W(y - 1)-^Z(G; (x - l)(y - 1), (y - 1)), (12) 

where we assume x ^ 1 and y ^ 1. 

We are going to use the following result on the complexity of exact evaluation of the Tutte 
polynomial. 

Theorem 4 (|23j). Exact evaluation of the Tutte polynomial is #P-hard for all rational numbers 
x, y except when 

1. (x - l)(y - 1) = 1; or 

2. (x,y) equals (1,1), (-1,-1), (0,-1) or (-1,0). 

The second part of Theorem [2] will be proved by reducing from exact evaluation of the Tutte 
polynomial. We prove the following Lemma in Section [3.11 

Lemma 3. Assuming the validity of the GRH, exact evaluation of the Tutte polynomial at x,y is 
polynomial-time Turing reducible to exact evaluation of R' 2 at A, fi, when 

(x- l)(y - 1) = 1/A- 1, y-l = Ai 2 , A {0, 1}, and fi^O. (13) 

Assuming GRH, by Lemma [3] and Theorem @J we have that exact evaluation of R' 2 at rational 
point (A, (j,) is #P-hard when A ^ {0, 1/2, 1} and \i ^ 0. We do not get #P-hardness for A = 1/2 
since the reduction is from the Tutte polynomial at (x — l)(y — 1) = 1 which is polynomial-time 
computable (part 1 of Theorem @|. (The other easy cases of the Tutte polynomial have no impact 
since y = 1 implies fi — and y G {0, —1} implies that /z is not real.) We proved the second part 
of Theorem [5J 

Now we prove the third part of Theorem [5] (the proof of main lemmas is deferred to later 
sections). To show #P-hardness of exact evaluation of R' 2 (G; 1/2, /j,) for [i ^ { — 1,0}, we prove a 
connection between R' 2 and the "permissive version of #BIS" (#PBIS) introduced in [17] : #PBIS 
is a generalization of #BIS where the weight of a set of vertices is determined by the number of 
pairs of neighboring vertices that are both in the set (in #BIS the weight is zero raised to the 
number of such pairs). 

^Permissive Bipartite Independent Sets with rj (#PBIS(r/)) 
Instance: a bipartite graph G = (U U W, E). 
Output: the quantity 

#PBIS(G;?7)= Yl {l + ri) w{rT) {l-rj) m ~ v>{ ' T \ 

<T:t/UlV->{0,l} 
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where w(a) is the number of edges in E with both endpoints labelled 1 by a. 

(We are using a different parametrization than [17] — our n and their 7 are connected by j 2 = 
(1 + 77)/(l - ry).) 
Note that 

#BIS(G) = 2 |B| #PBIS(G; -1). 

The following result is a generalization of Theorem[JJand shows that R' 2 encodes #PBIS(r/) as well. 
The proof is deferred to Section I3~2"1 

Theorem 5. Let G = (V, E) = (U U W, E) be a bipartite graph. 

#PBIS(G;r?) = 2^(G; 1/2, -r?). 

Note that #PBIS(G; 1) = 2l £ l +t , where t is the number of isolated vertices of G (since the other 
vertices have to be labeled 1 by a). This implies (TTTT) . 

The following result on #PBIS(?7) will be proved in Section [3731 

Lemma 4. #BIS is polynomial-time Turing reducible to #PBIS(?y) with r\ a rational number and 
77^ {±1,0}. 

The third part of Theorem [5] follows from Theorem [SJ Lemma 0] and the fact that exact com- 
putation of #BIS is #P-complete [53]. 

3.1 Reducing Tutte polynomial to R' 2 polynomial (proof of Lemma [3]) 

We will focus on bipartite graphs G — (U U W, E) such that vertices in partition W have degree 
at most 2 (a natural operation that produces such graphs is 2-stretch, that is, replacement of each 
edge with a path of length 2). 

Let G = (U U W, E) be a bipartite graph with max-degree in W bounded by 2. We call a 
connected component C = {Uq U Wc, Ec) of G pure if every vertex in Wc has degree 2 in C. A 
component that is not pure will be called mixed. The evaluation of R' 2 polynomial in G can be 
expressed using pure connected components as follows. 

Lemma 5. For every bipartite graph G = (U U W, E) such that the degree of each vertex in W is 
bounded by 2, 

fi;(G;A,,) = ^^'( s V S| , 

SCE 

where k'(S) is the number of pure connected components in (UUW,S). 

Before proving LemmaOwe need the following characterization of the rank of bipartite adjacency 
matrices over F2. 

Lemma 6. Let G = (U U W, E) be a connected bipartite graph with max-degree in W bounded by 
2. Let B be the adjacency matrix of G. Then 

, f \U\ if there is a vertex of deqree 1 in W , 

rk 2 (B) = \ ' , • 

I I (7 1 — 1 otherwise. 
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Proof. Let x G ¥ 2 be a solution of the linear system x 1 B — 0. Let Ui be the set of vertices u G U 
such that x u = i, for i = 0, 1. Note that no vertex v £ W has neighbors both in Uq and Ui 
(otherwise (x T B) v = 1). Thus for G to be connected either x = or x = 1. If there is a vertex 
of degree 1 in W then a; = 1 is not a solution and hence rk 2 (i?) = \U\. On the other hand if all 
vertices in W have degree 2 then x = 1 is a solution and hence rk 2 (L?) = |t/| — 1. □ 

Now we prove Lemma [5] 

Proof of Lemma\5i We will show that 

rk 2 (S) = \U\ - k'(S), (14) 

where rk 2 (5) is the rank (over F 2 ) of B, the bipartite adjacency matrix of (U U W, S), and k'(S) is 
the number of pure connected components of {U U W, S). 

Note that B has a block structure with a block for each connected component. The rank is the 
sum of the ranks of the blocks. Equation ([14")) now follows from Lemma [5J □ 

We now lay groundwork for the proof of Lemma [3] We use the following construction in the 
reduction. Given a graph H = (Vh,Eh) and a bipartite graph T = {Uy U W~r, E-f) with a specific 
vertex u £ Uy,we construct a bipartite graph G from H and T as follows. Let H = (U £ U W^, E^) 
be the 2-stretch of H, where corresponds to the vertices of H. For each vertex « in 0^- we 
identify v with u in a copy of T (thus we have \Vh\ = \Ug\ copies of T). We call the graph G the 
stretch-sum of H and (T,it). Note that if Wt contains only vertices of degree at most 2 then the 
partition of G containing Wjj contains only vertices of degree at most 2. 

We define two functions related to R' 2 . 

Definition 5. Let A, [i € K.. Let T ~ (U yjW, E) be a bipartite graph with a specific vertex u G U. 
Assume that the max-degree in W bounded by 2. We define 

Z' p (T;X^) = Y / ^ K ' iS) ^ Sl , (15) 

s 

where the sum is over all S C E such that u is in a pure connected component of (U U W, S), and 
k'(S) is the number of pure connected components of (U U W, S). 
Similarly, we define 

Z; n (T;A, M ) = ^A- K '( s y s l (16) 
s 

where the sum is over all sets S C E such that u is in a mixed connected component of (U U W, S). 

The following lemma provides a connection between the random cluster partition function Z of 
H and the R' 2 polynomial of G for rational A and [i. 

Lemma 7. Fix rational A {0, 1} and rational fi ^ 0. Let p be a prime such that A G Z*. Lei 
T = (J7 U W,E) be a bipartite graph with a specific vertex u G U, such that the max-degree in W 
bounded by 2. Suppose T satisfies 

X := XZ'JT; A,u) mod p, 

pV (17) 

Y := Z' m (T;\,fi) + XZ'JT; A,/i)=0 mod p. 
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Let G be the stretch-sum of H = (Vh,Eh) and (Y,w). Then 

R' 2 {G-, A, n) = \\ Vh \-\ u \x\ Vh \Z(H;1/\- 1,/i 2 ) mod p, (18) 
where Z(H; 1/A — 1,/i 2 ) is defined in JT]). 

Proof Let H = (VfrEg) be the 2-stretch of H. Note that H is a subgraph of G = (V G ,E G ). Let 
So Q E{j. Let As be a family of subsets S C .Eg such that S" n = So- Now, we evaluate 

2 A- K '( s y s i (19) 

seA So 

modulo p. 
Claim 1. 

__ _ J (A -1 — 1) K ( S °) [p s °>X' Vh ' mod p if (Vfj, So) has no mixed connected component, 
(|19p = < 

I mod p otherwise, 
where k(Sq) is the number of connected components of (Vfj, So)- 

Proof of Claim{^ Equation (|19[) can be rewritten as a product, where each term in the product 
corresponds to a connected component of (Va, So) (since each connected component with the copies 
of T attached to it influences k'(S) independently). Thus 

c 

where for each connected component C = (Vc,Ec) in (Vfr,So) such that there are k copies of T 
(we refer to the copies Ti, . . . , and to their special vertices u±, . . . , Uk) attached to it, 

$ c= ■■■ E A- K ' (U ^w s<UEc) /x |Bc|+E « W l^l. (20) 

SiCB Tl S fc CB Tfc 

Let Aifi be the set of Si such that m is in a mixed component of (VVi, Si) and let A^i = 2 Et * \ A^o- 
Equation ([20]) can be written as follows 



$ c =£;...£; J2 ••■ E A-*' (U «eM s<UBo V^ l+E<6wls<l - (21) 

xi=0 x k =0 S 1 £A 1 ^ 1 S k eA k ,m k 

We have 

«' ( |J SiUEc) = + 

\ie[fe] / ie[fc] »=i 

where I = 1 if a;i = • • • = Xk = 1 and C = (Vc,Ec) is a pure connected component of (V^,So), 
and I — otherwise. Thus 

52 52 A" K ' (U ^w 5 - u ^) M l £; cl+E ieW |s l | = 

SiGAi.a,! S k eA ktmit 

k (22) 

a-v |£cI E ••• E n A ~ K ' i(Si)+a; > |Sil = A-v |£cl ^ fe '(^-^)^', 

SieAi j!C1 S k eA k>Xk i=i 
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where k 1 = x% + • • • + Xk ■ 

Plugging ([22j into (JST) we obtain 



$ c = ^ Ec W k + L, (23) 

where L = (1/A— l)^ E °^X k if C is a pure component of (V H , So) and otherwise. Evaluating (j2"3"j) 
modulo p (using (|17p ) we obtain 

_ f mod p if C is a mixed component of (V h ,Sq), 

c = I (1/A - l)^ E °\X k modp otherwise. 

Thus (|19p is zero modulo p if there is a mixed component C in (V^pSo). Assume now that all 
components of (V^,So) are pure. The total number of edges in the components is |So|, the total 
number of copies of T in the components is \Vh\, and hence 

Ull) ee (1/A- 1) k ( s o) m IS | x |Vh mo d p . 

□ 

Now we use Claim Q] to prove (|18|) . Note that by Lemma [5j 

J^(G;A, M ) = Al^l 2 J2 A- K ' (S V |S| , (24) 

S C_E A SeAs 

where n'(S) is the number of pure connected components of (Vg, S). 

Note that by Claim[T]if So contains a mixed component then the inner sum in (|24j) is modulo p. 
Thus to evaluate (|24[) modulo p it is enough to sum over So which contain only pure components. 
Each such So is obtained from exactly one S' C Eu by 2-stretching. Note (V^,So) has the same 
number of connected components as (Vh, S'), and |So| = 2|S'|. Thus 

R' 2 (G;\fi) eee J2 - lr^'V 315 ' 1 

= X IVh1 ^X^ h1 Z(H;1/X- l,fi 2 ) mod p. 

□ 

We use different T for different values of fi in the reduction. When // ^ —2 we let Ti be the 
bipartite graph with bipartition U — {i*o, u±}, T4 7 = {uj | < i < k} and k + 2 edges: edge {uo, t>o}, 
and an edge between u\ and each Uj, for < i < fc. The specific vertex of Ti is uq. By elementary 
counting, we have 

AZ;(T i; A, /i) = (/i + l) fc+1 + /i 2 + A" 1 - 1, 
AZ;(T i; A, M ) + Z' m (T i; A,/x) - (At + l)((/x + l) fc+1 + A" 1 - 1). 

When /i = — 2 we let T2 be the bipartite graph with bipartition U = {uo, iti, W2}, W ~ {vi | < 
i < 2k} and 4k + 2 edges: {uq, vq}, {u±, «o}, and a complete bipartite graph between U \ {uq} and 
W \ {vo}. The specific vertex of T2 is uq. By elementary counting, we have 

XZ' (T 2 ; A, -2) = A~ 2 + S^A" 1 - 3 + 3 ■ 5 2k + A -1 , 
A^(T 2 ; A, -2) + Z^(T 2 ; A, -2) = -A~ 2 - S^A" 1 - 1 + 5 2k + 3A" 1 . 
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Fix rational A ^ {0, 1} and /i ^ 0. We want to find sufficiently many primes such that there is 
some integer k for which (f25j) satisfies (|17l) (when /.i ^ —2) or ([26]) satisfies p7|) (when /i = —2). 
We need the following result on the density of primes. 

Lemma 8 ([29, 30, 35 ). Let r,q G Q* and q ^ ±1. TTie density (inside the set of all primes) of 
primes p such that 

q k = r mod p, 

can be satisfied for some integer k is a positive constant, assuming GRH. 

Lemma [8] immediately yields the following two corollaries. 

Corollary 1. Fix rational A {0, 1} and rational fi {0, —2}. The density (inside the set of all 
primes) of the primes p such that there is an integer k for which (|25j) satisfies (|17|) is a positive 
constant, assuming GRH. 

Proof. If /j, = — 1, then to make (|2"5|) satisfy (|17p . it is sufficient to have A -1 ^ mod p. Thus, for 
all but a constant numbers of primes, and for all positive integers k, (l2~5j) satisfies (fTT|). 
Now assume fi ^ — 1. To make (|25p satisfy p7|) . it is sufficient to have 

(ii + l) k+1 + A" 1 - 1 = mod p, , 

' (27) 
\x ^ mod p. 

The corollary follows from Lemma [5] (and the fact that the number of primes such that /i = 
mod p is finite) . □ 

Corollary 2. Fix rational A g' {0,1} and rational /i = —2. TTie density (inside the set of all 
primes) of the primes p such that there is an integer k for which (j2l))) satisfies (fTTf is a positive 
constant, assuming GRH. 

Proof. We claim that if 

(A~ 2 - 3A" 1 + 1)(1 - A^ 1 )" 1 = 25 fe modp, 

, , (28) 
1 - A ^ 25 k mod p 

then (121)1) satisfies (|17l) . The first equation in (f2"5)l makes the second equation in (fT7| satisfied; and 
the first and the second equation in (|2"5)) make the first equation in (TT71) satisfied. 

By Lemma [8] the density of primes that make the first equation in (1281) satisfied is positive. 
Solving 

(A~ 2 - 3A" 1 + 1)(1 - A" 1 )" 1 = 1 - A" 1 mod p, 
we obtain A -1 = mod p and hence the second equation in (f28f is automatically satisfied. □ 

Proof of Lemma\^ Let x, y, A, \x be rational numbers such that (|13|) is satisfied. Suppose A = a/b 
and /i = c/d with a, b,c,d G Z*, gcd(a, 6) = 1, and gcd(c, d) = 1. 

Suppose we want to evaluate the Tutte polynomial for H = (Vh,Eh) at x, y. Let n :— \Vh\ 
and m := \Eh\- By (p~2|) . to evaluate T(H;x,y), we can instead evaluate Z(H; 1/A — 1,/U 2 ) (note 
that A 7^ 1 implies a; ^ 1 and y ^ 1 and hence (|12p applies). Recall that 

Z(H;l/X-l,^)= J2 (l/A-ir (S V |S| = ^r, (29) 
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where 

L= (b-ar^a n - K ^c 2 ^d 2m - 2 ^. 

SCE H 

Note that L 6 Z and |L| < 2 m |6 - a|"a"c 2m d 2m . 

We now prove the case fi ^ —2. For the case [i = —2, the proof is similar (by using T 2 and 
Corollary HJ. 

We choose n 3 primes pi , . . . , p n 3 such that 

• a, 6, c, d ^ mod ft, and a + b ^ mod ft, for each i 6 [n 3 ]; 

• there is some integer k for which (|27[) is satisfied with p = ft, for each i £ [n 3 ]; 

• ft = 0(n 4 ) for i g [n 3 ]; and 

• lYiliPi > 2 m+1 \b- a\ n a n c 2m d 2m . 

By Corollary[T] these primes exist. We can find them in time polynomial in n by exhaustive search. 

For each pi, let < h < pi be an integer for which ([25]) satisfies (fT7|) with p — pi (by Fermat's 
little Theorem, if ki is a solution, then ki + t(pi — 1) is a solution as well, for every t 6 Z). Again, 
we can find them in time polynomial in n by exhaustive search. 

We use Ti with k — ki as above. Let Gi be the stretch-sum of H and (Ti,uo) as above. Note 
that Gi has size polynomial in n since ki = 0(n 4 ). By Lemma[7]and (J29J), we have 

W = a n d 2m X- 2n X- n R^(Gi; A, /z) mod ft, 

where X = AZp(Ti;A,/i). We can make a query to the oracle to obtain the rational number 
R' 2 (Gi] A, fj.) and thus can compute L mod pi in polynomial time for each i € [n ]. By the Chinese 
remainder theorem, we can compute L in time polynomial in n (see, e.g., [3], p. 106). □ 

3.2 Reducing t^PBIS to the R' 2 polynomial (proof of Theorem [5]) 

Now we show the connection between #PBIS and the R' 2 polynomial; the proof of Theorem [5] is 
similar to the proof of the high-temperature expansion of the Ising model (see, e.g., |26j). 

Proof of Theorem \5\ 

#PBIS(G;t7) = (l + v) w{a) (l-v) lEl ~ w(a) 

ct:UUW^{0A} 

E II (l + VXHu),a(v))), (30) 

er:[/UVK->{0,l} {u,v}£E 

where 

v w v " I —1 otherwise. v ' 

Let 

*S,(ri,ff2 : = ]J X(c r i('«),c r 2(w)) and ^g,^ := E *s>i^2- 

{u,v}eS a 2 :W^{Q,l} 
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Expanding the product in (|30[) and changing the order of summation yields 
» = E E II VxHn),a(v)) 

rr:UUW^{0.1} SCE {u.v}GS 

= E^' E II X(cr(u),a(v)) 

SCE a:UUW^{0,l} {u,v}GS 

= E^' E (32) 

SCE CT1 :U-J-{0,1} 

Let Ns(v) denote the set of neighbors of v in the subgraph (UUW, S). Fix 5 and <j\ : U —> {0, 1}. 
We say that a pair S, <J\ is good if for every v £ W the number of vertices u £ N$(v) such that 
<7i (it) = 1 is even. A pair which is not good will be called bad. 

Claim 2. ^s,t?i = 2' H/ ' (— l)' s ' if the pair S, o\ is good; and ^s.o-i — if the pair S, o~\ is bad. 

Proof of Claim\^ Suppose that o 2 and a' 2 differ only in the value assigned to v £ W. For every 
u £ Ng(v) we have 



x(Mu),a 2 (v))-[ _ x(£7l(u)i ^ (u)) ifCTl(u) = 1 

Hence, 



* = f if |{u£iV s (u)|(Ti(u) = l}|isodd, 

*s>i,oi | * S , CT1 ,^ otherwise. 

If the pair S, cr\ is bad, then there is a vertex v € W such that \{u £ N$(v) | cti(m) = 1}| is odd. 
We can partition the W — > {0, 1} mappings into pairs <?2, o~' 2 that differ only in the label of v. For 
each pair, we have ^s,ai,<r 2 + ^S,ai,a' 2 = 0j an d thus $ s )CT1 = 0. 

If the pair S, o\ is good, then each ^s,(n,a 2 contributes the same value (— Since there are 
2^1 mappings from W to {0, 1} we have W 5j(Tl = 2^1 (-1) 151 . □ 

Claim 3. Fix SCE, the number of o\ : U — > {0, 1} such that the pair S,ai is good is 2l' 7 l -rk2 (' s ), 
where rk^S) is the rank (over ¥2) of the bipartite adjacent matrix of (U U W, S). 

Proof of Claim{^ Let B be the bipartite adjacency matrix of (U U W, S). Note that the pair S, o\ 
is good if and only if 

ajB = mod 2 

(we view o~\ as a vector with (a±) v = a\(v)). The claim follows from the fact that the number of 
vectors a £ {0, l}! 17 ! suc h that a T B = mod 2 is 2\ u ^ rk ^ B \ □ 

By Claim [2] and Claim H (j32j) equals 

J2 2 iq+|w|-ri«(S) ( _^|5| = 2 \ U \+\ W \R' 2 (G; 1/2, -rj). 

SCE 

□ 
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3.3 Reducing #BIS to #PBIS (proof of Lemma SJ) 

We use the following construction in the proof of Lemma |U Given a graph G = (V,E), a prime 
p and a positive integer k we construct a bipartite graph G' by replacing each vertex v G V by a 
cloud of fcp new vertices and each edge e G £ by a cloud of p — 1 new vertices. For every vertex 
v G V and every edge e adjacent to v we add a complete bipartite graph between the cloud of v 
and the cloud of e. 

Lemma 9. Let p > 2 be a prime, k be a positive integer, and rj be a rational number such that 
r/ G Z* and 

1 + 77 



+ 1 = mod p, 
l-r) J (33) 

— 1 ^ mod p. 

Let G = (V,E) be a graph and let G' — (V',E') be the graph constructed above (using p,k). Then 
the number of independent sets of G is congruent to $:PBIS(G"; rj) mod p. 

Proof. We are going to evaluate 

#PBIS(G';r ? )= E J] (1 + VX (a(u),a(v))) (34) 

er:V'^{0,l} {ti,»}EE' 

modulo p (where x is given by (|3"Tj) ). 

Let s be a vertex of G. Consider a set of configurations (that is, er's in ([34]) ) that differ 
only in the labels assigned to the vertices in the cloud of s (thus the configurations agree in the 
labels assigned to the vertices outside of the cloud of s). Note that by symmetry the value of 
ri{u v}eE'(l + ? 7X( cr ( M )j f7 ( w ))) depends only on h, the number of vertices in the cloud of s that are 
assigned 1. There are ( k £) such configurations and hence the total contribution of configurations 
with h G {1, . . . , kp — 1} to (|3"4")l is zero modulo p. Thus, we only need to consider the cr's such that 
for every s € V the vertices in the cloud of s are assigned the same label. 

For q : V — ► {0, 1} let C ? be the set of assignments V — > {0, 1} that assign label ?(s) to every 
vertex in the cloud of s (for every s G V) . Now we evaluate 

e n a +»«(*(«), <*>))). (35) 

ueC s {»,«}6E' 

Note that each vertex in the edge clouds influences its own set of terms in the product in (|35|) . 
Thus the sum in (|35[) turns into the following product 

63= II *«s),<;m (36) 

where 

y) = [(1 + TO (x, 0)) fc P(l + r? X (y, 0)) fcp + (1 + VX (x, l)) kp (l + vx(y, l)) fcp ] P_1 . (37) 
Now we evaluate \l/(x, y) for cc, y G {0, 1} 
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• Case 1: x = y = 1. We have 

*(1,1) = [(l-r,f kp + (l + v) 2kp ] P ~\ 
and from the assumption (|33|) we obtain 



#(1,1) = (1 - r,) 2kpi - p -V 



2kp 



p-1 



— mod p. 



(38) 



• Case 2: x = y = 0. We have 

#(0,0) = [2(1 -r)) 2kp ] P ~ 1 = 1 modp, (39) 
where in the last congruence wc used r\ ^ 1 mod p and Fermat's little Theorem. 

• Case 3: x = 1, y = (the case a; = 0, y = 1 is the same). We have 

*(1,0) - [(1 -7/) 2fcP + (l+ 77)^(1 -7 7 ) fcp ] P_1 



(1 - ■q) 2kp ^ p ~ 1 ^ 
1 mod p, 



I + 77 



fcp 



p-i 



(40) 



where the last congruence follows from Fermat's little Theorem and the fact that ( ( 1 + 77) /(l - 
rf)) k ^ —1 mod p (since otherwise ((1 + rj)/(l ~ rj)) 2k = 1 mod p which violates ([33])). 



From (|38|) . (|39|) . and (|40j) we obtain that (|36|) modulo p is 1 if c corresponds to an independent 
set of G and is zero otherwise. Thus we have 

#PBIS(G';t7)= II *(?(«),?(*)) = #BIS(G0 mod p. 

?ev^{o,i} {s,t}e£ 

□ 

Fix rational 77 ^ {±1,0}. We want to find sufficiently many primes p such that (|33|) can be 
satisfied for some integer k. The following results shows that the density (inside the set of all 
primes) of such primes is a positive constant. 

Lemma 10 ([!]). For q € Q* and q =/= ±1, the density (inside the set of all primes) of primes for 
which there exists fc S Z such that 

q k = — 1 mod p, 

is a positive constant. 

The mapping r\ \-t (jz^j) 2 maps ryeQ*\{±l}toQ*\{±l} and hence we have the following. 

Corollary 3. Fix rational rj {±1,0}. The density (inside the set of all primes) of the primes p 
such that (|33p has an integer solution to k is a positive constant. 
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Proof of Lemma^ The proof is a routine application of the Chinese Remainder Theorem (similar 
to the proof of Lemma [3]). Given an instance G — (V,E) of #BIS, let n := \V\ and L be the 
number of independent sets of G. Note that L < 2™. 
We choose n 2 primes p%, . . . , p n i > 2 such that 

• for each i £ [n ], (j53")l has an integer solution to k with p = pi; 

• Pi = 0(n 3 ) for i £ [n 2 ], and 

By Corollary [3j these primes exist when 77 ^ {±1,0}, and we can find them in time polynomial in 
n using exhaustive search. 

For each p if let < k, t < p { be an integer solution of ([3"3")l with p — Pi (if h is a solution of ([551) 
then fcj +i(pi — 1) is a solution for every t £ Z), again we can find the fcj using exhaustive search. We 
then construct a bipartite graph G\ as above from G, pi and fc^. Note that G\ has size polynomial 
in n since ki < pi — 0(n 3 ). As (l33|) is satisfied, by Lemma [9j we have 

L = #PBLS(<^; rj) mod p t . 

We can make a query to the oracle to obtain the rational number #PBIS(G^; rj) and thus can 
compute L mod pi in polynomial time for each i £ [n 2 ]. By the Chinese remainder theorem, we 
can compute L in time polynomial in n. □ 

4 The linear-width of a graph 

For the proof of Theorem[3jwe will use the linear-width of a graph, a concept which was first defined 
by Thomas |37j . In this section we prove a bound on linear- width in terms of tree- width. 

The linear-width of a graph G = (V, E) is the smallest integer £ such that the edges of G can 
be arranged in a linear order ei, . . . , e m in such a way that, for every i £ [m], there are at most 
£ vertices that have an adjacent edge in {ei, . . . , e^-i} and an adjacent edge in {e^, . . . , e rn }. It is 
known that computing the linear- width of a graph is NP-complete [3B]. For paths and cycles the 
linear-width is easy to compute. 

Example 1. The linear-width of a path is 1. The linear-width of a cycle is 2. 

Let ei, . . . , e m be a permutation of the edges of G — (V, E). We say that a vertex v £ V is 
dangerous w.r.t. i £ [m], if there exist two edges e^e^ adjacent to v such that j < i < k. Let Di 
be the set of vertices which are dangerous w.r.t. i. Note that the linear-width of G is the minimum 
value of maxi \Di\ optimized over all permutations of the edges. 

Now we give an upper bound on the linear- width for trees. 

Lemma 11. Let T = (V,E) be a tree on n vertices. The linear-width of T is at most \}og 2 n\ . 

Proof. We describe an ordering a witnessing that the linear- width is at most Ll°g2 n \ ■ Pick an 
arbitrary vertex r £ V to be the root. For a vertex v, let T v be the subtree rooted at v and let 
\T V \ be the number of vertices in T v . We do a depth-first search (DFS) on T where the children are 
explored in the increasing order of \T W \ (that is, the smaller subtrees are explored first). Let a be 
the order in which the edges are discovered by the DFS. 
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Next, we show that the linear- width of a — ei,...,e m is at most [log 2 tt-J ■ Let i £ [to]. Let 
vi, . . . ,vj be the vertices on the stack right before the edge e,{ was discovered (thus is adjacent 
to Vj). Note that the vertices which are not on the stack are not dangerous w.r.t. i (they either are 
already fully explored or haven't been discovered yet). Note that if Vk+i is the rightmost child of 
Vk then all the edges adjacent to Vk are discovered and hence Vk is not dangerous w.r.t. i (for all 
k 6 \j — 1]). On the other hand if Vk+\ is not the rightmost child of Vk then our exploration order 
implies 

\T Vk+1 \ < \T Vk \/2. (41) 

Note the following two facts: Vj may be dangerous and \T V . | > 2. If there are £ dangerous vertices 
among vi, . . . , Vj-i then by |4T|) we have \T W \ > 2 • 2 l where w is the topmost dangerous vertex. 
Using \T W \ < n we obtain the result. □ 

For general graphs we will show a generalization of Lemma 111! a bound on the linear- width of 
G in terms of the tree- width of G. We now define tree- width (see, e. g., [27] for a nice treatment). 

Given a graph G = (V,E),& tree decomposition of G is a pair (T, {Uh}hev T ) where T = (Vt, Et) 
is a tree and Uh Q V satisfy: (i) each edge of G is in at least one subgraph induced by Uh', and 
(ii) for any three vertices t\,t 2 , t 3 of T such that t 2 is in the path between ti and t$ in T we have 
U tl H U t3 C U t2 . The width of a decomposition is max^evr Wh\ ~ L The tree-width of G (denoted 
tw(G)) is the minimum width optimized over all tree decompositions. 

Lemma 12. Let G = (V, E) be a graph. Then 

linear-width(G) < (tw(G) + l)(|_log 2 n| + 1). 

Proof. Given a graph G with \V\ = n, let (T, {Uh}hev T ) (where T = (Vr, Et)) be an optimal tree 
decomposition of G. Let r :— \Vr\- We can assume that r < n (see, e. g., [37] p. 579). Let a' be the 
order in which the vertices of T are discovered by the DFS of Lemma [TT1 W.l.o.g. let a' = 1, . . . , r. 

We will need the following "separation" property of dangerous vertices: a path between an 
explored and an unexplored vertex must contain a dangerous vertex (this follows from the fact that 
an explored vertex with an unexplored neighbor is dangerous). 

Now we construct the ordering a = e\, . . . , e m on the edges of G. For each h £ Vr, let Gh be 
the induced subgraph of G on Uh- An edge e € E may appear in several Gh] we "assign" it to the 
first Gh (in the order a'). Let a be an ordering such that edges assigned to Gh come before the 
edges assigned to Gh 1 for any h < h! . Note that we do not impose any particular order on the edges 
assigned to Gh (for any h G [r]). 

Let i € [to] and let h be such that ei is assigned to Gh- Let Di C V be the set of dangerous 
vertices w.r.t. i. Let D' h C Vt be the set of dangerous vertices just after h was discovered. Let 
v E Di\Uh- We will show that v £ U g for some g £ D' h . 

Let Cjj and e^ 2 be adjacent to v and such that i\ < i < Note that and e^ 2 cannot be 
assigned to Gh (since v is not in Uh). Thus e+ x is assigned to an explored vertex ji of T and ej 2 
is assigned to an unexplored vertex j 2 of T. Note that v is in Uj 1 and Uj 2 and hence it is in Uk 
for all k that are on the path between j\ and j 2 in T. By the "separation" property of dangerous 
vertices one of them (say g) must be dangerous, and by the second property of tree decomposition 
U g contains v. 

Hence, 

\D i \<\U h \ + ]T \U g \ < (tw(G) + l)(Llog 2 nJ+l). 

□ 
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5 Analysis of the single bond flip chain for trees 

Given a tree G — (V,E), let be the set of 2^1 subsets of E. By Lemma [2 for every H C E, 
we know that rk 2 (H) is the size of maximum matching of the subgraph (V, H). Let w(H) be the 
size of maximum matching in a graph (V,H). Let P be the transition matrix of the single bond 
flip Markov chain M. from definition [2] It's easy to see that Ai is ergodic with unique stationary 
distribution n such that 

n(H) oc A™( H V |B| - 
The goal of this section is to prove Theorem [3] 

5.1 The canonical paths 

We will bound the mixing time of our chain Ai using the canonical paths method, introduced in [8] 
[34l [25]. Now we go over the basic definitions for Markov chains, see, e. g., [24] for a comprehensive 
background. 

Definition 6. The total variation distance of two probability distribution v and v 1 on is 

Hen 

Definition 7. The mixing time from initial state H , th{s), is defined as 

r H {e) = min{t : \\P\H, ■) - n\\ TV < e}, 
and the mixing time r(e) of the chain is defined as 

r(e) = max{r H (£)}. 

/it" 

Let a — e\, . . . , e m be an ordering of the edges of G — (V, E) (we will usually use the orderings 
supplied by Lemma [TT] or Lemma [T2"]) . Given any pair /, F £ O, let I © F be the symmetric 
difference of / and F (that is, the set of edges which are in either I or F but not in both). We 
define a canonical path 77^ between I and F as follows. Let , . . . , ei k be the edges from I © F 
ordered according to a (that is, i\ < %% < • ■ ■ < ik)- Let 

1ItF = {H 0> H x ,...,H k ), (42) 

where H = I, H k = F and Hj = Hj-i © {e^}. 

Lemma 13. Let G = (V, E) be a graph. Let a = e%, . . . , e m be an ordering on E with linear-width 
I. Let L,F be subsets of E and let H be on the canonical path (|42l) (that is, H = Hj for some 
j G {0,...,k}). Then 

\w(I) + w(F) ~ w{H) - w(C)\ < £, 

where C = L © F © H . 

Proof. Let Q = {ex, . . . , e^}. Note that H = (F n Q) U (L n Q c ), where Q c is the complement of Q 
(that is, E\Q). Similarly, C = (L n Q) U (F n Q c ). 
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Let D be the set of dangerous vertices w.r.t. +1. Let Mj and Mp be maximum matchings of 
I and F, respectively. Let 

M H = {M F n Q) U {Mi n Q c ) and M c = (Mj n Q) U (Mp n Q c ). 

Note that all vertices of Mjy with degree > 2 are in 73 (a vertex which is not D has all adjacent 
edges (in G) from Q or from Q c and hence the adjacent edges (in Afe) agree with Mj or Mj?). The 
same is true for Mc- Moreover if a vertex v £ D has degree 2 in M# then it has degree in Mc- 
Thus by removing < \D\ edges from Mh and Mc we can turn both of them into matchings. Thus 

w(H) + w(C) > w(I) + w(F) - \D\ > w(I) + w(F) - i. (43) 

Note that a canonical path from V := H to F' := C passes through H' := I (with C := I' ® F' © 
H' = F). Thus 

«;(/) + w(F) = w(H') + w(C) > w(l') + w(F') -£= w(H) + w{C) - t. (44) 
Combining (l43l) and p4|) we get the lemma. □ 



5.2 The congestion of M. 

Now we analyze the congestion of the collection L = {ji,f \I,F £ £1} where 'ji.f are canonical 
paths defined in (g2J|. For each transition (H, H') such that P{H. H') > 0, let cp(H, H') be the set 
of pairs (I, F) such that (H, H') £ Ji,p- The congestion of F on (if, iJ') is (see, e. g., [24] ) 

= pJhJf) 2^ -^HT bl ' Fl (45) 

V ; I,F:(H,H')e-yi,F 

where |7/,f| is the length of Ji,f- The congestion of L is defined as 

g:= max Q(h,h>)- 

P(-ff,/J"')>° 

We will use the following connection between the congestion and the mixing time. 

Theorem 6 (0134]). t h (s) < g(log(l/7r(ff)) + log(l/e)) for each starting state H 6 0. 

At the end of this section we prove the following bound on the congestion of T. 

Lemma 14. Let G = (V, E) be a graph. Let a = e\, . . . , e m be an ordering on E with linear-width £. 
For every (H,H') such that P(H 1 H') > 0, and for every A,/i > we have 

q {h , w) < 2\E\ 2 \ e , 

where A = max{A, 1/A}. 

We can now prove Theorem [3] 

Proof of Theorem^ Since G — (V, E) is a tree, by Lemma [Til we have I < [log 2 , by Lemma [T4l 
we have 

Q < 2\E\ 2 \ e < 2|£:| 2 n |lo fe A| < 2n 2+|log ^ A| . 
Theorem [3] now follows from Theorem [6] □ 
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Now we bound the congestion of our canonical paths. 

Proof ofLemma[TJ\ We will bound q {h ,h>) for every (H,H') such that P(H,H') > 0. Let H = H 
if n(H) < tt(H') and H = H' otherwise. Note that 

^ = n(H)P(H, H>) = ir(H')P(H', H), (46) 

since M. is reversible. We define a mapping / : cp(H, H') — » such that /(/, F) — I F H for 
every pair (I,F) £ cp(H,H'). 

First, note that / is an injection. Given J £ Q we can determine the unique /, F such that 
/(/, F) = J, by first computing J © H , and the using the ordering a on the edges of G to recover 
/ and F. 

Note that 

|/| + \F\ = \H\ + \f(I,F)\, (47) 

and 

\w(I) + w(F) - w(H) - w(f(I, F))\ < I, (48) 

where (|48|) follows from Lemma fl~3l 

Let L = Y.j \ w{J) ^ A - We have the following upper bound on q^h.h 1 )- By (|45jl and (|46|) . we 
have 

(I,F)ecp(H.H') 7r ^ n ) 



X w(I)+w(F)-w(H)^\I\ + \F\-\H\ 

L 



2\E\ 2 E 

(FF)£cp(H,H') 



< 2|F| 2 ^ E * ^ (49) 

(I,F)£cp(H,H') 

< 2\E\ 2 X i , (50) 
where (|4"9")l follows from (|4"7T) and (|4"5)l . and (I5TH) follows from the fact that / is an injection from 

c P {h, h') to a □ 



6 The single bond flip chain for the random cluster model 

The following sampling problem corresponds to the random cluster model defined in |T|). 

Random Cluster Model with q,n>Q (RC(g,/i)) 
Instance: a graph G = (V, E). 
Output: S CE with probability of S oc q K ^^. 

The random cluster model was introduced by Fortuin and Kasteleyn [15] as a model in statistical 
physics, and was intensively studied since then (see, e.g., the monograph by Grimmett [21] )■ By 
the results of [231 S21 SI] ) the problem of exactly computing the partition function of RC(q,/z)) is 
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T^P-hard even in bipartite planar graphs for all rational numbers q, fi > except when (i) q = 1; 
and (ii) q = 2. 

In the context of approximate counting, little is known about ^RC(q, fi). The case q — 2 and 
H > corresponds to the ferromagnetic Ising model and has an FPRAS [26]. For "dense graphs" 
an FPRAS exists for all q, fi > 0, [2 CQ. 

Cooper and Frieze [7] and Cooper et al. [B] investigated the Swendsen-Wang algorithm for 
RC(g,/x). They proved, using coupling and conductance methods, that the Swendsen-Wang algo- 
rithm is rapidly mixing on some classes of graphs. Gore and Jerrum [20 pointed out that the single 
bond flip chain for random cluster model is torpid (mixes in exponential time) on the complete 
graph K n for q > 3 and a value of /i depending on the size of the graph (that is, /x is not a fixed 
constant). 

We will show that the single bond flip chain mixes for graphs with bounded tree-width. Our 
motivation is "showing that Markov chain mixes" , not approximation of #RC((j, /i) (since the prob- 
lem of computing the partition function of RC(g,/i) is polynomial-time solvable |31j for bounded 
tree- width graphs). 

The single bond flip for the random cluster model is analogous to the chain from Definition [5J 
At time t, pick an edge e £ E at random and let S = X t © {e}. Set X t +i = S with probability 
(1/2) min{l, g K ( 5 )~ K ( x *)^l s H^*l} ancl X t +i = X t with the remaining probability. The chain M 
defined above is ergodic with unique stationary distribution n such that tt(H) oc q K ( H ^ ij) h \ . Let P 
be the transition matrix of A4. 

We will bound the mixing time of M. for every fixed constant q, fi > 0. Given an ordering a 
on the edges of G, we define the same canonical paths for each pair of subgraphs /, F, and the 
same mapping / for every (H, H') such that P(H, H') > 0, as in Section [5j We present a result 
analogous to Lemma IT3l 

Lemma 15. Let G — (V, E) be a graph. Let a = e±, . . . , e m be an ordering on E with linear-width 
I. Let I,F be subgraphs of E and let H be on the canonical path (|42j) (that is, H = Hj for some 



where C = I © F © H. 

Proof. As in the proof of Lemma [TBI let Q = {e±, . . . ,ej.}. Let D be the set of dangerous vertices 
w.r.t. ei +i. We can split V \ D into V\ U V2, where V\ contains vertices whose adjacent edges are 
from Q and Vi contains vertices whose adjacent edges are from Q c . 

Let Ii be the subgraph of / induced by Vi U D and Fi be the subgraph of F induced by Vi U D 
(for i = 1,2). Note that H can be viewed as a gluing of I2 and F\ on D; we will denote this 
H = I 2 UF 1 . Similarly C = h U F 2 . 

First note that 



since F± , F 2 can by obtained from F by splitting each vertex in D and the splitting of one vertex 
increases the number of components by at most one. 

Let /' consist of the connected components of I that do not contain a vertex from D. Let I[ be 
the subgraph of / induced by Vi (for i — 1,2). Note that 



j G {0,...,A;};. Then 



+ k(F) - k(H) 



k(F 1 )+ K (F 2 )<k(F) + \D\, 



k(Fx U I' 2 ) + k(F 2 U /() = + k{I' 2 ) + k(F 2 ) + = 

k(Fi) + k(F 2 ) + k(I') < k(F) + \D\ + k(I') < k(F) + k(I) + \D\, 



(51) 
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where the first and second inequalities follow from the fact that the following pairs of graphs are 
vertex disjoint: F\ and L 2 , F 2 and I[, and I[ and I' 2 . 

Note that I2 © I 2 consists of connected components that contain a vertex from D. Thus k(F\ U 
h) = k(F% U 1 2) (since adding a component from h © I' 2 to Fi U I' 2 only attaches that component 
to an existing component). Similarly k(F2 U ii) = r(-F2 U /{). Combining with (|51j) we obtain 

«(#) + k(C) < k(J) + k(F) + \D\ < k(I) + k{F) + t. (52) 

Note that a canonical path from I' := H to F' := C passes through H' := I (with C := I' © F' © 
H' = F). Thus 

k(7) + k(F) = K(fr') + k(C) < «(/') + k(F') +£ = k{H) + k(C) + 1. (53) 

Combining (f52|) and (|53|) we get the lemma. □ 

We have the following bound on the congestion. 

Lemma 16. Let G = (V, E) be a graph. Let a = e±, . . . , e m be an ordering on E with linear-width 
t. For every (H,H') such that P(H,H') > 0, and for every q,/j, > we have 

Q(H.H') < 2\E\ 2 q e , 

where q = max{q, 1/q}. 

Proof. Same as the proof of Lemma [14J □ 

By Theorem [6] and Lemma 1161 we have the following result. 

Theorem 7. Given an graph G = {V, E) with m — \E\, n = \V\ and linear-width £, the mixing 
time r(e) of A4 is 

r(e) = O (m 2 q 4 (m + n\ logg| + m\ log^| + log(l/e))) , 

where q = max{g, 1/q}. 

By Lemma [T2l we conclude the following. 

Corollary 4. Let k be an integer and let q, ji > 0. Let G — (V, E) be a graph with tree-width 
bounded by k, n := \V\, and m := \V\. The mixing time t(e) of the single bond flip chain M. for 
HC(q,fi) is bounded as follows 

r(e) = O (m 2 g (fc+1)(1+lo ^ n) (m + n\ logg| + m| log/*| + log(l/e))) , 

where q = max{g, 1/q}. 

We note that we do not need to find the linear- width (or tree- width) of a graph efficiently, the 
existence result is enough for our proof of the mixing time. 
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7 Conclusions 



We conclude with an observation that a generalization of RWM(A,/z) does not have an FPRAS 
(unless RP=NP) and a few questions. 

Let A be an m x n matrix whose entries are zeros, ones, and indeterminates, where each 
indeterminate occurs once. A completion of A is a substitution of 0, 1 to all the indeterminates in 
A. We denote Ca to be the set of all completions of A. Let rk2(B) be the rank of B over F2. Can 
we sample B from Ca with the probability of B proportional to A rk2 ^ ? Note that this problem is a 
generalization of the RWM(A, 1) problem. It turns out that finding the minimum rank completion 
of a matrix is NP-hard (Proposition 2.1, [32]) and hence a sampler is unlikely (unless NP=RP), 
since for A = 2 _Tl a random completion will be the minimum rank completion (with constant 
probability). The sampling problem could be easy for sufficiently large A (the problem of finding 
maximum rank completion is in P, see, e. g., Section 4.1 of [22]). 

Question 1. What other interesting properties are encoded by the polynomial? 

Question 2. Can one sample maximum rank completions of a matrix? 
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